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FOREWORD 


The  influence  of  a  background  plasma  on  the  diocotron  instability  of  a 
relativistic  hollow  electron  beam  is  investigated  within  the  framework  of  a 
cold  fluid  model.  The  background  plasma  is  taken  to  have  a  uniform  density 
except  at  the  mean  radius  of  the  electron  beam  where  the  density  is  assumed  to 
be  discontinuous.  Of  the  charged  species  comprising  the  plasma,  only  the  electrons 
influence  the  dynamics  of  the  beam-plasma  system;  the  effect  of  the  ions  is 
negligible  because  of  their  small  plasma  frequency.  A  dispersion  relation  for 
the  eigenfrequencies  of  the  system  is  derived  and  employed  to  examine  instabilities. 
It  is  found  that  the  fundamental  (Z  =  1)  mode,  which  is  unconditionally  stable  for 
a  hollow  electron  beam  in  a  vacuum,  is  now  unstable  for  certain  parameter  regimes. 
The  appearance  or  growth  of  instabilities  with  increases  in  the  plasma  density 
discontinuity  is  observed,  while  increases  in  the  relativistic  factor  y  are  seen 
to  exert  a  stabilizing  influence.  It  is  also  found  that  thin  beams  are  generally 
more  unstable  than  thick  ones. 
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INTRODUCTION 


Numerous  devices  employing  relativistic  hollow  electron  beams, 

1  2  3 

e.g.,  the  coaxial  autoaccelerator,  the  free  electron  laser,  ’ 

4  5-7 

the  relativistic  magnetron,  and  the  gyrotron  are  currently  objects 

O 

of  study  in  diverse  areas  of  research.  A  recent  investigation0 
examined  the  stability  of  relativistic  hollow  electron  beams  in  a 
vacuum.  In  this  paper  the  investigation  of  stability  is  extended  to 
include  the  case  in  which  a  background  plasma  is  present.  This  case 
is  of  interest  because  any  operating  experimental  device  such  as  those 
listed  above  will  contain  a  plasma  produced  by  the  ionization  of  the 
residual  gas  by  the  electron  beam.  If  the  density  of  this  plasma 
is  sufficiently  great,  its  presence  may  have  important  consequences 
for  the  operation  of  the  device  since  the  beam  parameters  required  for 
stable  operation  may  be  affected. 

The  analysis  of  the  beam-plasma  configuration  will  be  carried  out 
within  the  framework  of  a  macroscopic  cold  fluid  model,  the  temperatures 
of  the  beam  electrons,  plasma  electrons,  and  plasma  ions  all  being 
assumed  negligible.  Computation  reveals  that  the  ions,  by  virtue  of  their 
relatively  small  plasma  frequency,  can  be  neglected  in  the  perturbation 
analysis  employed  in  examining  the  stability  of  the  beam-plasma 
configuration;  the  ions  do  affect  the  equilibrium  state,  however. 

The  analysis  is  restricted  to  the  case  where  the  electron  beam  is 
tenuous  and  where  the  electron  density  of  the  background  plasma  is  no 
greater  than  that  of  the  electron  beam.  In  specifying  a  form  for  the 
density  of  the  background  plasma  It  is  presumed  that  in  a  physical  device 
this  density  would  exhibit  a  gradient  across  the  interior  of  the  electron 
beam.  For  the  sake  of  mathematical  tractability  this  gradient  will  be 
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accounted  for  by  assuming  that  the  plasma  density  is  sharply  discontinuous 
at  a  region  inside  the  electron  beam;  elsewhere  the  plasma  density  is 
assumed  to  be  uniform. 

We  shall  begin  by  determining  the  equilibrium  state  of  the 
beam-plasma  configuration.  The  stability  properties  of  the  system 
will  then  be  investigated  through  the  use  of  a  perturbation  method 
employing  the  Maxwell-fluid  equations.  In  particular  we  shall 
derive  a  dispersion  relation  for  the  eigenfrequencies  of  the  perturbed 
configuration.  Next  this  dispersion  relation  will  be  employed  to 
investigate  the  stability  properties  of  the  configuration  for  various 
ranges  of  the  parameters.  It  will  be  seen  that  the  fundamental  ( i  *  1) 
mode  can  be  unstable,  in  contradistinction  to  the  vacuum  case  where  this 
mode  is  unconditionally  stable.  It  will  also  be  observed  that 
increases  in  the  plasma  density  discontinuity  or  decreases  in  the  beam 
thickness  are  generally  destabilizing,  whereas  increases  in  the  rela¬ 
tivistic  factor  Yb  have  a  stabilizing  effect. 
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BEAM-PLASMA  CONFIGURATION'  AT  EQUILIBRIUM 

Figure  1  illustrates  the  configuration  of  the  electron  beam  and 
the  background  plasma  at  equilibrium.  With  reference  to  the  cylindrical 
coordinate  system  shown  in  Figure  1,  the  beam  is  seen  to  be  propagating 
with  uniform  velocity  B^ce^  and  is  coaxial  with  a  perfectly  conducting 
outer  cylinder  of  radius  R^.  R^  and  R2  are,  respectively,  the  inner  and 
outer  radii  of  the  hollow  electron  beam.  A  uniform  axial  applied 
magnetic  field,  Bgez>  is  assumed  to  be  present  to  provide  radial 
confinement  of  the  beam. 

The  equilibrium  electron  density  profile  of  the  beam  and  of  the 
plasma,  nj^(r)  and  n^(r)  respectively,  are  illustrated  in  Figure  2. 

Both  profiles  are  assumed  to  be  independent  of  0.  The  beam  electron 
density  is  taken  to  be  uniform  across  the  width  of  the  beam  and  to 
fall  sharply  to  zero  at  the  beam  edges.  The  plasma  electron  density  is 
taken  to  be  piecewise  uniform,  undergoing  a  discontinuous  decrease  by 
a  factor  a  at  the  midpoint  of  the  electron  beam  from  its  value 
at  the  beam  axis.  The  plasma  is  assumed  to  be  electrically  neutral 
at  equilibrium  so  that  each  ion  density  profile  of  the  various  ion 
species  would  have  the  same  form  as  that  of  the  plasma  electrons. 

The  equilibrium  electron  densities  may  be  written  as 

n°(r)  =  nbg[U(r  -  R^)  -  U(r  -  R2) ]  ,  (1) 

n®(r)  =  neg[b'(r)  -  (1  -  a)U(r  -  Rp)]  ,  (2) 

where  R  *  (R.  +  R«)/2,  where  U(r)  is  a  unit  step  function  (Heaviside 

p  1  2 

function),  where  n^g  end  ngg  are  constants,  and  where  the  parameter  a 
satisfies  0  <  a  <  1. 
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In  analyzing  the  equilibrium  configuration  the  following  additional 
assumptions  are  made: 

i)  All  of  the  beam  and  plasma  variables  are  azimuthally  symmetric, 
i.e.,  3/30  =  0.  The  azimuthal  motion  of  the  beam  will  thus  be 
laminar. 

ii)  The  beam-plasma  configuration  is  infinitely  long  with  no 
axial  variations,  i.e.,  3/3z  =  0. 

iii)  The  beam  and  plasma  charged  particles  are  in  slow  rotational 

x  £  drift  velocities. 

iv)  The  azimuthal  momenta  of  the  charged  particles  are  non- 
relativistic  and  small  compared  to  the  axial  momentum  of  the  beam  electrons. 

v)  The  beam  electron  density  is  assumed  to  be  such  that  the  following 
inequalities  are  satisfied: 


equilibrium  with  mean  azimuthal  velocities  equal  to  £ 


(3) 


/<1)  .<*>  <<  U  , 

cj  ce  pb 


(4) 


2  2 

Here  w  ,  =  4ire  n,  »/v,  m 
pb  bO  'b  e 

frequency,  where  =  1 
cyclotron  frequency;  <uce  -  Ybucb  t*ie  P^asma  electron  cyclotron 
frequency;  and  <ocj  i  q^Bg/mjC  is  the  cyclotron  frequency  of  the  jth 
species  of  ions,  where  qj  is  the  ionic  charge.  Inequality  (3)  implies 
that  the  electron  beam  is  tenuous,  while  (4)  states  that  the  beam  is 
nevertheless  sufficiently  dense  so  that  the  lower  hybrid  frequency  of  the 
plasma  is  much  less  than  the  plasma  frequency  of  the  beam  electrons. 

As  stated  earlier,  the  analysis  will  be  carried  out  within  the 
framework  of  a  cold  fluid  model.  The  equations  of  motion  and  continuity 
which  cl  -acteriz*"  he  Pth  component  of  the  fluid  (beam  electrons, 


is  the  square  of  the  beam  electron  plasma 

/A  -  B?;  to  ,  =  eBn/y,m  c  is  the  beam  electron 
b  CD  u  d  e 


10 


NSWC  TR  81-393 


plasma  electrons,  or  the  jth  species  of  ions)  are  then 


+  V 

'Viy 


V)y  m  V 
P  P^P 


ewVl  +  C%V  *  V 


(5) 


3_ 

3t 


n 

P 


+  7 


(n  V  ) 
p^p 


0  , 


(6) 


where  £  and  respectively,  the  electric  and  magnetic  fields;  where 

n^  and  ^  are  the  number  density  and  velocity,  respectively,  of  the 

fluid  component;  where  and  are,  respectively,  the  mass  and 

charge  magnitude  of  the  particles  comprising  the  fluid  component; 

where  is  the  sign  of  the  charge  of  the  fluid  component  particles; 

2  -1/2 

and  where  y  is  the  relativistic  factor  (1  -  V  •  V  /c  )  for  the 

y  'V-y  *Vy 

particles  of  the  fluid  component  (y  is  essentially  unity  except  for  the 
beam  electron  component) . 

The  electromagnetic  field  is  governed  by  Maxwell's  equations 


7  *  Jg 


—  —  B 
c  3t  ^ 


(7a) 


V  x  B 

'V 


+  —  -I—  E 

c  3t  ^ 


(7b) 


7  •  E  =  4ttp  ,  (7c) 

where  p  and  J  are  the  charge  and  current  densities,  respectively. 

The  electric  and  magnetic  fields  can  be  derived  from  potentials  <j>  and  A, 

*  ■  -V*  -“i  .  (8) 

<» 
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These  potentials  satisfy 


=  -4irp  , 


Ha 

h  1  H j . 
\,  =2  3t2r  c  ' 


(10a) 


(10b) 


At  equilibrium,  consistent  with  the  assumptions  stated  earlier, 
we  have 


0. 

n  =  n  (r)  , 
b  b 


(11) 


V 

'V/p 


v  „(r)&  +  g  ce  , 
p0  'VQ  p  -vz 


(12) 


=  4>  (r)  , 


(13) 


£  *  Az(r>!2  ’ 


(U) 


where  the  superscript  zero  denotes  the  equilibrium  form  of  the  associated 
quantity.  The  following  also  apply  at  equilibrium: 


p  *  -en 


b  ' 


(15) 


Z  *  "6bcenb£z  * 


(16) 


0  . 


where  n^  is  given  by  Eq.  (1).  By  employing  the  technique  of  Uhm  and 

g 

Siambis  the  following  expressions  are  obtained  for  the  equilibrium 

rotational  velocities  oi®  =  V°_/r  of  the  electrons  in  slow  rotation: 

b  b® 

u)£(r)  =  uD[l  -  (R1/r)2]U(r  -  , 


(17a) 


for  the  beam  electrons,  and 
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wp(r)  s  ^buD{U  '  <Vr)2lU<r  "  Rl) 

-  [1  -  (R2/r)2]U(r  -  R2)>  , 

for  Che  plasma  electrons,  where  Is  Che  diocotron  frequency, 

WD  5  “pb/2Vcb  * 


(17b) 


(18) 


It  can  also  be  shown  that 


u>?  -  0(u°)  ,  (19) 

J  P 

for  the  plasma  Ions  in  slow  rotation;  the  explicit  form  of  the  ion 
rotational  velocities  will  not  be  needed. 
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The  stability  of  the  electron  beam-plasma  configuration  will  be 
examined  by  applying  a  perturbation  analysis  to  linearized  equations  of 
motion  and  continuity.  In  the  perturbation  analysis  a  typical  quantity 
of  interest,  e.g.,  charge  density  p,  is  written  as 

p(r,t)  =  p°(r)  +  6p(r,t)  ,  (20) 

where  p^  is  the  equilibrium  charge  density  and  Sp  the  perturbed  charge 
density.  In  the  stability  analysis  all  of  the  perturbed  quantities 
are  Fourier  decomposed  in  a  manner  typified  by 

6p(£,t)  *  6pA(r,6)exp[i(Jl0  -  wt)]  ,  (21) 

where  l  is  the  azimuthal  harmonic  number  and  ui  is  the  eigenfrequency . 

Using  the  perturbed  quantities  Sp^,  etc.,  a  dispersion  relation  will 
be  derived  for  the  eigenf requencies .  By  solving  this  dispersion  relation 
for  the  eigenfrequencies  the  stability  properties  of  the  configuration 
can  be  determined. 

The  first  step  in  the  perturbation  analysis  is  the  determination  of 
the  perturbed  electric  and  magnetic  fields  6E,  6B.  These  fields  can 
be  derived  from  perturbed  potentials  6$,  5A  in  the  manner  shown  in  Eqs.  (8) 
and  (9)  (except  that  the  constant  term  B^e^  in  Eq.  (9)  is  now  dropped). 

The  perturbed  potentials,  together  with  perturbed  source  terms,  satisfy 
relationships  given  by  Eqs.  (10a, b). 

The  stability  analysis  is  restricted  to  frequencies  satisfying 

|«oR  |  <<  t.c  .  (22) 

With  this  restriction  the  following  are  obtained  from  Eqs.  (10a, b) 


15 


NSWC  TR  81-393 


for  Che  perturbed  potentials: 


(k  i_  r  L.  _  i! ) 

lr  3r  r  3r  ^2 J 

(l  3_  3_ 

l  r  3r  r  3r  “  r2y 


6^i(r,w)  *  -^itSp^r.w)  , 


5Azl(r,oi)  -  -  —  <5Jzi(r,oj) 


(23a) 


(23b) 


In  Eq.  (31b)  and  <SA z  are  the  z-components  of  the  perturbed 

current  density  and  vector  potential,  respectively.  Because  the 
azimuthal  velocities  of  the  components  of  Che  beam-plasma  configuration 
are  assumed  to  be  nonrelativistic,  the  forces  resulting  from  the 

9 

magnetic  field  generated  by  the  azimuthal  current  are  negligible. 

Hence  only  the  axial  components  of  the  perturbed  current 
density  and  perturbed  vector  potential  are  needed  in  the  analysis. 

By  virtue  of  the  foregoing  considerations  the  perturbed  electro¬ 
magnetic  field  can  be  obtained  from  the  perturbed  potentials  as 
follows: 


<5Era(r,w)  = 


“574*i(r’“)  » 


6E0a(r,a») 


IV 

-  —  (r ,<d)  , 


(24b) 


5BrJl(r,u>)  -  —  5AzJl(r,w)  , 


(24c) 


6B9*(r’“) 


r—  6A  .  (r ,u>) 
3r  z  i  ’ 


The  equations  of  motion  and  continuity  for  the  perturbed  fluid 
variables  are  obtained  by  linearizing  Eqs.  (5)  and  (6)  about  the 
equilibrium  values  of  these  variables,  neglecting  terms  higher  than 
first  order.  The  results  for  the  equation  of  motion  are 
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-i(oj  -  iuP)6V  -  (2oj^  +  oj  )6V.  . 

y  rux  y  cy  6y£, 


- ^  I"  -  8  6A  )  , 

Y  m  3r  i  \izl 

u  u 


(25a) 


I  u°  +  | —  (roj^)  +  oj  ]5V  .  -  i(co 

I  u  3r  V  cyj  rpi 


(25b) 


-  lu  )<5V 

u  8pJl 


—  (6$  -  g  6A  .  , 

Y  m  r  y£  %  zir  * 
u  u 


while  for  the  equation  of  continuity  the  following  is  obtained: 


fl  Of i  a 

-i(w  -  i.w  )<5n  +  n  —  —  (rSV  .) 

\i  V  U(r  3r  ryJt 


(26) 


+  —  5V 
r  8p£ 


+  (|-  n°)6V  „  =  0 

3r  u  ryi 


In  these  expressions  the  subscript  u  indicates  the  component  of  the 
configuration  (beam  electrons,  plasma  electrons,  or  one  of  the  ion 
species);  the  symbol  denotes  the  cyclotron  frequency  of  the  component. 

The  perturbation  solution  proceeds  by  solving  Eqs.  (25a, b)  and 
(26)  for  the  perturbed  number  densities  in  terms  of  the  perturbed 
potentials.  By  employing  the  resulting  expressions  for  the  perturbed 
number  densities  in  the  relations 


5p, 


e  q  5n 
UP  V 


(27a) 


6Jzi  *  -eSbc<5nb  ’  (27b) 

we  can  obtain  expressions  for  the  perturbed  charge  and  current  densities 
in  terms  of  the  perturbed  potentials. 
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The  expressions  for  the  perturbed  charge  and  current  densities 

in  terms  of  the  perturbed  potentials  are  employed  in  Eqs.  (23a, b) 

in  order  to  obtain  a  pair  of  coupled  differential  equations  for  the 

perturbed  potentials.  These  differential  equations  can  be 

considerably  simplified  by  employing  inequalities  (3)  and  (4).  In 

particular,  the  effects  of  the  ions  can  be  shown  to  be  negligible 

in  comparison  to  those  of  the  beam  and  plasma  electrons.  The  influence 

of  these  latter  two  components  depends  principally  on  terms  with  the 
3  2  9  2 

coefficients  —  u>  ,  ,  —  w  where 
3r  pb  dr  pe 

“pb  S  4,te2nb/Ybme  ’  (28a) 

w2  =  4ire2n^/m  ;  (28b) 

pe  e  e 

use  of  Eqs.  (1)  and  (2)  then  leads  to 

h  “Jt  ■  v[5<r  ■  v  • 6(c  *  V1  •  <29*> 

fe  v  • -  V  •  <29b> 

«  2 

where  S>^  is  the  square  of  the  plasma  frequency 

ii2  =  4ire2n  n/m  .  (30) 

pe  eO  e 

After  defining  the  generalized  potentials 

«*,  *  «♦,  -  «b«*„  .  <»*> 

“t 5  8b5*i  -  vA«  •  «lb> 

the  following  set  of  coupled  differential  equations  are  obtained 
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DISPERSION  RELATION 


A  dispersion  relation  for  the  eigenf requencies  of  the  beatn-plasma 
configuration  can  be  derived  by  solving  Eqs.  (32a, b)  after  imposing 
the  boundary  conditions  that  <5^,  65  are  finite  at  r  *  0  and  vanish  at 
r  *  R^.  The  derivation  of  this  relation  will  now  be  outlined  briefly. 

First,  we  note  that  the  right-hand  sides  of  Eqs.  (32a, b)  vanish 
except  at  r  3  R^,  R^,  R2;  for  r  not  equal  to  any  of  these  values  the 
solutions  for  6^,  6?^  will  both  have  the  form  arm  +  a'r  m  where 
a, a*  are  independent  of  r.  Therefore,  the  r-interval  (0,Rc>  is  divided 
into  four  subintervals  (O.R^),  (R^,  R^) ,  (Rp,  R2),  (R2>  Rc> ,  and  the 
solutions  for  6!^,  65  in  each  of  these  subintervals  are  written  as 

*  akr*  +  a^r“*  (33a) 


6? 


I 


+ 


(33b) 


where  the  subscript  k  indicates  the  subinterval  (k  *  1,2,3, 4). 

Now  the  potentials  6?^  are  continuous  functions  of  r  at  the 

points  r  ■  R^,  Rp,  R2-  The  r-derivatives  of  the  potentials  are  dis¬ 
continuous  at  these  points  by  quantities  which  are  obtained  by  integrating 
Eqs.  (32a, b)  with  respect  to  r  over  infinitesimal  intervals  centered 
at  these  points.  For  instance,  the  integration  of  Eq.  (32a) 
from  r  *  Rj  to  r  *  r|  yields 


3__ 

3r 


6i^(r,u)) 


(r,w) 


2  luir 


RjU  -  a‘%(R1) ] 


5 (R1 ’ * 

04) 


By  using  this  and  similar  results  at  r  »  R^,  Rj,  and  by  employing  the 
boundary  conditions,  we  obtain,  after  some  tedious  algebra,  the  dispersion 
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relation 


d3(ai/wD)3  +  d2(u/a)D)2  +  d1(u/u)D)  +  dQ  -  0  ,  (35) 


where 


d3  ~  ferf-Da-V  > 


(36a) 


d2  3  fS12  '  A(5c2  *  1)(X2  +  VJ(1  ~ 


cp 


(36b) 


+  -  »  +  <1  -  V<c2  '  V1  ■ 


dl  =  [l  (5c2  “  1)ApA2  ~  4?12Ap  "  *(S12  "  5c2)A2 


Cc2512H1  ‘  +  YbAa[BbecpU12 


-  tuc2  -  1)X7}  +  -  uc 


(36c) 


c2<’P2  "  ‘^c2  *  ’p2)A2 


-  (5c2-S2)ClpKl-  V]  * 

d0  3  l^h2  *  5c2)A2  +  ?C2  5!23(1-  «cp>Ap 

+  Yb'at,£(«c2  '  Ep2»2  •  5c2Ep2Hl 


-  Ecp>Elp  +  8bV‘<Ee2  -  E12>»2 


5c2512^  ’ 


(36d) 


and  where 
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!w»  !  1  -  ‘VV"  •  (37,) 

X2  E  wb(R2)/wD  ’  <37b) 

»p  *  “p«p)/"D  • 

Xa  =  (l-oJOpo/i^Q  .  (3?d) 

It  is  significant  to  note  that  the  plasma  parameter  a  and  the 
equilibrium  electron  densities  enter  into  the  dispersion  relation 
coefficients  solely  through  the  parameter  X^  given  by  Eq.  (37d). 

The  other  parameters  in  these  coefficients  depend  only  on  the  geometric 
quantities  R^,  R2,  Rp,  Rc>  and  on  the  relativistic  quantities  yb» 
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RESULTS 


The  dispersion  relation  given  by  Eq.  (35)  has  been  employed  to 
study  the  stability  of  the  electron  beam-plasma  configuration  for  a 
substantial  range  of  parameter  variation.  In  the  following  are 
presented  the  most  significant  results  from  this  study. 

One  of  the  first  observations  was  the  appearance  of  an  instability 
in  the  l  »  1  mode  for  certain  parameter  regimes.  This  contrasts 
sharply  with  the  unconditional  stability  of  this  mode  in  the  absence 

Q 

of  plasma.  In  Fig.  3  is  plotted  the  region  of  instability  in  the 

(R./R,,,  X  )  parameter  space  for  y.  -  2  and  R./R  =  0.8.  The 
i.  t  a  b  L  c 

growth  in  the  area  of  the  unstable  range  of  with  increasing  X^  is 

evident.  It  is  also  observed  from  additional  numerical  studies  that 
increases  in  y^  cause  the  unstable  region  to  shrink,  and  that  generally, 
the  thinner  the  electron  beam,  i.e.,  the  closer  R2/R2  is  to  unity, 
the  larger  the  unstable  range  in  Xa-  This  last  property  is  found  to 
be  true  for  all  modes. 

It  is  also  of  interest  to  examine  the  stability  in  the  (R^/R^,  R2/Rc> 
parameter  space  for  a  fixed  X^  for  the  l  -  1  mode.  Numerical  studies 
show  that  the  appearance  of  an  instability  is  usually  more  sensitive 
to  changes  in  the  ratio  R^/R2  than  to  changes  in  the  rat:i-0- 

It  is  again  found  that  the  unstable  region  shrinks  with  increases  in  y^. 

The  decreases  in  the  areas  of  the  unstable  regions  in  the  (R. /R0,  x  ) 
and  (R^/Rj,  Rj/I^)  parameter  spaces  with  increases  in  y^  appear  to  be 
unique  to  the  l  -  1  mode.  For  the  higher  modes  the  converse  behavior 
is  generally  the  rule. 

The  shape  of  the  unstable  region  in  the  (R^/^,  R2/Rc)  parameter 
space  for  the  £  -  2  mode  is  found  to  differ  considerably  from  those 
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of  the  l  =  1  and  l  >  2  modes.  The  shape  of  the  unstable  region  for  i  =  2 

is  shown  in  Fig.  4;  the  relevant  parameters  are  y.  =  5  and  X  =0.5. 

o  a 

In  Fig.  4  it  can  be  seen  that  the  appearance  of  an  instability  can  be  very 

sensitive  to  the  R«/R  ratio. 

i  c 

The  areas  of  the  unstable  regions  are  found  to  decrease  as  the 

mode  index  increases.  In  Fig.  5  is  shown  the  shape  of  the  unstable 

region  in  the  (R-j/R2»  X  )  parameter  space  for  the  i  =  3  mode  and 

for  *  5,  “  0-8-  As  the  mode  index  is  increased  while  the 

other  parameters  (y^.  R2/Rc)  are  held  constant,  the  extreme  left-hand 

boundary  of  the  unstable  region  is  found  to  move  to  the  right. 

This  section  on  results  is  concluded  with  an  examination  of  the 

behavior  of  the  eigenfrequencies  with  parameter  variations.  In  Fig.  6 

are  plotted  the  real  and  imaginary  parts  of  the  normalized  (with  respect 

to  alp)  eigenfrequencies  for  several  modes  for  varying  X^.  The  parameters 

for  these  plots  are  *  5,  R^/R2  *  0.75,  and  R2/Rc  *  0.8.  It  is 

evident  that  the  growth  rates  of  the  instabilities,  i.e.  ,  the  imaginary 

components  of  the  eigenfrequencies,  increase  as  X  increases  for  all 

a 

except  the  i  ■  2  mode.  Thus  we  see  that  the  greater  the  plasma  electron 
density  discontinuity,  the  more  unstable  the  beam-plasma  configuration. 

The  effect  of  relative  beam  thickness  on  the  eigenfrequencies  is 
typified  by  the  plots  of  Fig.  7,  which  illustrate  the  variations  in 
the  growth  rates  of  the  instabilities  with  relative  beam  thickness  for 
beam-plasma  parameters  of  y^  =  5,  R2/Rc  a  0.8,  and  X^  =  1.0.  It  is 
evident  that  as  the  relative  beam  thickness  decreases  from  maximum 
(S, /R2  *  0)  the  number  of  unstable  modes  and  the  peak  growth  rates 
increase,  reach  maximum  values,  and  then  fall  to  zero  as  the  beam 
thickness  approaches  zero.  Thin  beams  are  seen  to  be  more  unstable. 
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in  general,  than  thick  beams. 

Finally,  additional  numerical  studies  show  that  while  normalized 

growth  rates  tend  to  increase  as  y^  increases,  the  absolute  growth  rates 

decline  by  virtue  of  the  normalization  factor  being  inversely 
2 

proportional  to  y^- 
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The  stability  of  a  relativistic  hollow  electron  beam  propagating 
through  a  background  plasma  with  a  step  discontinuity  has  been 
analyzed  within  the  framework  of  a  cold  fluid  model.  The  analysis 
began  with  the  solving  for  the  equilibrium  state  of  the  beam-plasma 
configuration  (Sec.  II),  after  which  a  perturbation  technique  was 
employed  with  linearized  equations  of  motion  to  derive  differential 
equations  for  generalized  field  potentials  (Sec.  III).  A  dispersion 
relation  for  the  eigenf requencies  of  the  configuration  was  then 
obtained  (Sec.  IV)  and  was  used  to  examine  the  stability  properties. 

The  significant  findings  include  the  appearance  of  an  instability  in 
the  l  ■  1  mode,  which  is  stable  in  the  absence  of  plasma;  the  generally 
more  unstable  nature  of  thin  beams  as  opposed  to  thick  ones;  the 
stabilizing  influence  of  large  values  of  -yb;  and  the  increase  in  the 
growth  rates  of  the  instabilities  (except  for  the  l  -  2  mode)  with 
increases  in  the  plasma  density  discontinuity. 
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CONDUCTING  CYLINDER 


Electron  beam-plasma  configuration  and  coordinate  system 
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ELECTRON  BEAM 
NUMBER  DENSITY  n° 

/ 

-  Rp  =  1/2  (R-,  +  R2) 


PLASMA  ELECTRON 
NUMBER  DENSITY  n! 


i 

i 


Electron  number  density  profiles  of  electron  beam  and  plasma 
at  equilibrium 
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Fig.  3  Instability  region  in  the  (R^/R^ ,  Xo)  parameter  space  for  the 
l  *  1  mode  for  Rj/Rg  "0.8  and  *  2 
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0.0  0.5  1.0 

Rn/R2 

Fig.  5  Instability  region  in  the  (R^R^  Xq)  parameter  space  for  the 
l  ■  3  mode  for  ®2^Rc  "  and  ■  5 
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Fig.  7  Dependence  of  the  normalized  growth  rates  on  relative  beam 

thickness  for  several  modes  for  y.  ■  5,  X  -  1,  and  R-/R  •  0 

D  L  C 
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